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^ ■ Abstract 

On 
O 

l/*-) ' method is shown as well as a role of approximate resonances in wave 

turbulence theory. Applications of clipping method are sketched for 
gravity-capillary and drift waves. Brief discussion of possible transi- 
£pH| tion from continuous spectrum (= kinetic equation) to discrete spec- 

^ ■ trum (= clipping) is given at the end. 

B 

> 

^ ■ Wave kinetic equation has been developed in 60-th (see, for instance, 

[7j, ^J) an d applied for many different types evolution PDEs. Kinetic 
equation is approximately equivalent to the initial nonlinear PDE but 
has more simple form allowing direct numerical computations of each 
wave amplitudes in a given domain of wave spectrum. Wave kinetic 
equation is an averaged equation imposed on a certain set of correla- 
tion functions and it is in fact one limiting case of the quantum Bose- 
Einstein equation while the Boltzman kinetic equation is its other 
limit. Some statistical assumptions have been used in order to obtain 
kinetic equations and limit of its applicability then is a very com- 
plicated problem which should be solved separately for each specific 
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equation [T7] . 



Simply formulated, conditions used for obtaining of kinetic equa- 
tion can be described as following: 

• each wave takes part in resonant interactions, 

• each wave interacts with all other waves simultaneously, 

• all wave amplitudes are of the same order. 

As a result, original evolution equation has been reduced to an 
equation of a form 



which is solved in respect to each separate wave amplitude A{. Here 
function Q depends on the form of initial nonlinear PDE and notation 
8(^2 ki) i s use d for delta-function which is equal to zero on the solu- 
tions of the system of equations describing exact resonances among 
interacting waves 



and is non-zero otherwise. Wave kinetic equation is often written 
out in the form 



because amplitude square is proportional to the wave energy and this 
form allows to treat the results produced by kinetic equation in terms 
of energy exchange between the interacting waves. 

Existence of resonances, i.e. solutions of Sys.|(U), in a (n + l)-wave 
system for some specific PDE allows us to replace the PDE by kinetic 
equation which govern energy transfer through the spectrum. A very 
interesting fact is that even some classification of dispersive PDEs due 
to their integrability properties was constructed basing on solutions 
of Sys.Q. We present it below very briefly and for more details see 



First of all, it was proven that in case when Eq.(^l) does not have 
any solutions, initial nonlinear PDE could be transformed into some 
linear PDE by canonic transformations. If, on the other hand, Eq.(^l) 





(1) 
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does have some solutions, then original PDE still has nonlinearity after 
canonic transformation and it can be written out as 

T i 8(k l ±k 2 ±...±k t ) 

i ~ ~ ~ — 

uj(ki) ± oo{k2) ± ... ± uj(ki) 

where 5 is a delta-function. Vertex coefficient Tj for some specific 
io is just generalization of interaction coefficient from the system for 
slowly changing amplitudes in case of io-waves interactions. In case 
of zero vertex coefficient, Tj = we fall into the class of equations 
having solitonic solutions and in case of non-zero vertex coefficients, 
as it was shown before, the kinetic equation is constructed which to 
some extent is equivalent to the original PDE. 

Kinetic equations have been written for infinite interaction do- 
mains, i.e. for an infinite plane or an infinite channel, and were used 
successfully for about 25 years for description of many types of waves, 
mostly in cases for 3- and 4-waves interactions. The results of lab- 
oratory experiments showed that "reasonable agreement is normally 
obtained with theoretical results for the infinite case" [2] , if the wave- 
length are small enough in comparison to the size of the experimental 
basin. The cases when wavelengths are comparable with the charac- 
teristic sizes of the experimental basin remain, as a rule, unexplained 
6 and have been named "effects of finite lengths" .2]. Attempts to put 
some additional physically relevant terms into the kinetic equations 
in order to make them applicable to long-wavelength systems have 
failed. The author of the pioneering work in this field, O. Phillips, 
who obtained the first kinetic equation in wrote in |12j that "new 
physics, new mathematics and new intuition is required" in order to 
understand energetic behavior of large-scale systems). 

Roughly speaking, large-scale wave system is a system where wave 
do notice the boundaries. It means that original PDE has to be re- 
garded with zero or periodical boundary conditions and correspond- 
ingly wave spectrum is discrete. Now Sys.Q is a system of Diophan- 
tine or algebraic equations to be solved in integers and assumptions 
used for obtaining of kinetic equations in case of continuous wave 
spectrum together with discrete Sys.Q have to be scrutinized closely. 
Why do any pair of integers (m, n) has to be part of some solution? 
Or, even more, to be part of infinitely many solutions? Is it really 
true that this algebraic system always has infinitely many solutions? 
Rather not. At least, one has to prove it. 

It was shown [B] , [HJ that the study of discrete systems demands a 
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new approach, different from method of kinetic equation, due to their 
quite different properties: 

• all interacting waves are partitioned into disjoint classes, there is 
no energy flow between these classes and it is possible to study 
each class independently; the number of elements of a class is 
not large, mostly classes consist of a few waves only; 

• there are many waves which do not take part in any resonant 
interactions at all (for instance, in case of atmospheric planetary 
waves this amounts to 60 % of all waves); 

• wave interaction is local in a sense that for any fixed wave its 
interaction domain (i.e. the one which contains all the waves 
which can interact with a given one) is finite and can be written 
out explicitly; 

• number of interacting waves depend drastically on the shape of 
a basin (for example, on the ratio of its sides); there exist many 
basin where wave interactions are not possible; for some fixed 
dispersion it is possible to describe set of all such basins. 

In particular, it means that discrete wave systems can not be de- 
scribed in terms of kinetic equation and has to be studied separately. 

All the results |Hj were obtained for exact resonances while in ^01 
it was shown that specific properties of discrete systems are still valid 
for some non-zero resonance width 

< n = u(ki) ±u{k 2 ) ± ... ±u(ki). (2) 

The existence of allowed resonance width gave rise in ^01 to formula- 
tion of Clipping method and these results we briefly present in Section 
2, just for completeness of exposition, which we briefly present in Sec- 
tion 2. Applications of Clipping method for two physical problems are 
demonstrated in Sections 3 and Section 4. Brief discussion is given in 
Section 5. 

2 Approximate resonances 

In order to answer the questions about resonance width, we need some 
estimation of as a linear form on the value of the function to taken 
in different points of its definition domains depends on the to changing 
domain. Let us consider cases. 

1. Let u : Z x Z — > Q, where Z denotes the set of integers and 
Q denotes the set of rational numbers. Then obviously f2 can be 
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represented as difference of two rational numbers a/b and c/d and we 
have trivial estimation 

l n l = l?-3l*H >0 ' (3) 

As illustrative example for this case one may take spherical Rossby 
waves with ui — 



n(n+l) ' 

2. Let w:ZxZ-» Q(o), where a is an algebraic number, i.e. it 
is a zero of some polynomial 

P(x) = aQX r + aix r_1 + ... + a r , 

where not all equal to zero. The field Q(ce) denotes the algebraic 
expansion of Q, i.e., the set of numbers of the form 

a c 
b + d a ' 

with a, b,c,d € Q. In this case we may use generalization of the 
Thue-Siegel-Roth theorem If the algebraic numbers a±, 02, a s 
are linearly independent with 1 over Q, then for any e > we have 

+ <?2«2 + ••• + q s a s ~p\> cq~ s ~ c , 

for all p, qi, q2, q s £ ^ with q = maxj \qi\. The constant c has to 
be constructed for every specific set of algebraic numbers separately. 
As illustrative example for this case one may take capillary waves with 
lo = (m 2 + ra 2 ) 3 / 2 . 

3. Let u:Zx2->l, where M denotes the set of real numbers, u is 
an arbitrary real-value function of integer variables. To find frequency 
discrepancy 0, in the finite spectral domain D with some T such that 

D = {(m, n) : < m, n < T < 00}, 
it is enough to calculate it as 



17 = min p Vtp, 

where 

n p = w(A*) ± ± ... ± 

fcj = (m^,np, Vj = l,2, ...s, 
w(j^) ± uj{k p 2 ) ± ... ± Lj(k p s ) + Vp, 
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and p is finite because the total number of wave vectors belonging 
to D is finite. So denned £l p obviously is a non-zero number as a 
minimum of finite number of non-zero numbers. As illustrative exam- 
ple for this case one may take gravity water waves with lo = k tanh ak. 

Therefore, for a wide class of dispersion functions u it is possible 
to find the low boundary for the frequency discrepancy fl; in a finite 
spectral domain this boundary exists for arbitrary u. 

Of course, in order to apply the theoretical results of wave tur- 
bulence theory for some specific physical problem, one has also to 
compute values of wave amplitudes corresponding to the regime of 
weak turbulence. In other words, one has to characterize explicitly 
the nonlinearity of a system under consideration, i.e. to choose a 
small parameter < e « 1 describing a weakly nonlinear regime. 
For instance, for planetary waves the ratio of the particle velocity to 
the phase velocity is often taken as the small parameter which makes 
it possible to obtain ^H] following estimation for a wave amplitude 
a(m,n) with a wave vector k = (m,n): 

6mn\2 2n+1 - m 

\Qj\Tn Ti)\ <C 

1 v ' M 5n(n+l) m + n+3 (n-m)(5n-m-3)' 

Usually it is a non-trivial task to compute small parameter e symbol- 
ically, and also technical facilities has to be taken into account (see 
Section 3) in order to get some plausible results when planning some 
laboratory experiment. 

It is important to realize that in order to observe effects predicted 
by wave turbulence theory, one has to compute together 1) allowed 
resonance width, and 2) allowed magnitudes of resonantly interacting 
waves. 



3 Gravity-capillary waves 

Now that we know exact results and are able to calculate some small 
but non-zero resonance width let us plan some laboratory experiment 
- an experimental procedure which allows laboratory experiments to 
be performed in order to observe how our theory works in reality. 
Water waves are probably the most appropriate class of wave systems 
for the experimental performance because their properties have been 
studied profoundly and the experimental facilities developed are very 
sophisticated indeed. We regard here 
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dispersion relation for gravity-capillary waves 



u\k)=gk + ^k 3 



v 



where h is undisturbed depth, g is the gravitational acceleration, 
v is the density, fi is the surface tension, k is the wave vector, k = \k\ 
and in c.g.s. units g=981, fi = 74 (for water "filtered of particles 
with nominal sizes larger than 0.2 \xm and with temperature about 
18° C", see 0, p. 55), v = 1. The ratio \ijv may vary considerably 
even for the same liquid, for instance, if the water is not filtered of 
particles then ^i/v = 54; this ratio may also depend on the tempera- 
ture [Sj • Therefore, first of all we have to choose the spectral domain 
of interest and then find resonantly interacting waves corresponding 
to the given ratio \ijv. Let us regard spectral domain k x ,k y < 30 
which corresponds to wavelengths in the range (0.1, 3.0 cm) and find 
corresponding solutions of resonance conditions in the form: 



Choice of frequencies to generate 

Let us notice that when theoretically we have to find only exact 
solutions, in practice all the solutions with small enough frequency 
discrepancy have also to be considered. The ratio D of frequency 
discrepancy to the minimal driving frequency can be chosen as the 
characteristic of the discrepancy smallness and the value of D is obvi- 
ously connected with experimental precision. Driving frequency can 
be controlled with precision of order 10 -5 |3]; therefore we may only 
consider as resonant triads with small enough D, say D = 10 -6 . Be- 
low we present a few solutions (three couples of numbers in square 
brackets denote the wave numbers of three wave vectors while three 
numbers in round brackets mean the three corresponding frequencies 
in Hz) for different liquids: 



TypeA 

li/v = 7h : [1,2] [9,1] [10, 3]; (8.7638, 40.4435, 49.2073) 

fj,/u = 47 : [1,26] [16, 4] [17, 30]; (147.0295,75.8317,222.8612) 

y,/u = 27 : [1, 10] [28, 6] [29, 16]; (30.7235, 129.5023, 160.2258) 

H/v = 16 : [1,6] [4, 5] [5, 11]; (15.5681, 16.2945, 31.8626) 
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We present here just a few solutions from the multitude of existing 
solutions in order to show that for wide variety of liquids resonant 
triads do exists (the value of \ijv = 75 corresponds to clear water, 8° 
C; n/v = 47: glycerine C 3 H 5 (OH) 3 , 20° C; p/v = 27: benzol C 6 H 6 , 
60° C; n/v = 16: benzaldehyde C 6 H 5 CHO film on water, 20° C). 

In order to demonstrate one of the most striking properties of dis- 
crete resonant systems - namely, the existence of many non-interacting 
waves - we have also to compute the frequencies of the the waves pro- 
viding big frequency discrepancy. To find them, it is enough to solve 
the system of equations 



D = max k ^ 3 D 123 , 

D 123 x min(u}(ki),u(k2),u{fa)) = 
= w(fci) +uj{k 2 ) -uj(k 3 ), 
kxi ~\~ k X 2 = k X j, , ky 1 + ky 2 = ky 3 , 
k Xi i ky i £ N, i = 1,2,3, k Xi , ky i < 30. 

As the solutions of this system we obtain the triads of waves with 
the maximal possible ration of frequency discrepancy to the mini- 
mal driving frequency in the given spectral domain k Xi , k yi < 30 and 
for given \ijv (it is enough indeed to find triads with discrepancies 
larger than experimental precision, i.e. to replace its first equation by 
D > 0.1. Below a few triads with large discrepancies are presented: 



TypeB 
(i/v = 75 
fi/v = 47 
lijv = 27 
fi/u = 16 



[11, 15] [14, 15] [25,30]; (112.6460, 130.0788,337.7987) 
[14, 14] [15, 16] [29, 30]; (98.6504, 114.4728,295.8396) 
[4, 4] [26, 26] [30, 30]; (16.2595, 186.8502, 230.8321) 
[5, 5] [25, 25] [30, 30]; (17.8606, 137.0759, 178.8991) 



Thus now we have two types of triads: type A (exact resonant 
triads) and type B (triads with big discrepancies). Excitation of the 
frequencies corresponding to A-triads provides the possibility to ob- 
serve standard periodic energy exchange within the triad. Excitation 
of the frequencies corresponding to some .B-triad will not provide any 
periodic motion for similar initial conditions (i.e. for the same exper- 
imental facilities, the same liquid and the same magnitudes of initial 
wave amplitudes). 

Choice of initial amplitudes 
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Now we have found the frequencies of waves to generate but we 
still have to choose initial values for the wave amplitudes. The reso- 
nant interaction theory deals with weak nonlinearity, i.e. some small 
parameter < e << has to be chosen in order to estimate how big 
the initial amplitudes are allowed to be: they have to be big enough in 
order to eliminate the linear wave propagation but at the same time 
the amplitudes have to be not too large in order to escape turbulence. 

There exist different ways of choosing e, the choice is defined by 
the specifics of the wave system. For instance, as it was mentioned 
above, for spherical planetary waves e is chosen as a ratio of the parti- 
cle velocity to the phase velocity while for water waves usually e = ak 
, where a is amplitude of a wave and k = \k\ . This quantity character- 
izes the steepness of the waves and the values e = 0.1 or 0.2 correspond 
to the weakly nonlinear regime. Therefore, for arbitrary given wave 
numbers it is easy to calculate the appropriate wave amplitudes. 

Dependence on basin form 

Now we are going to behavior difference in triad's behavior in dif- 
ferent experimental basins. Suppose we a rectangular experimental 
basin; each specific mode has the form 

r a ,m\ „ /mX mirx niry 

[A mn (T) cos — + B mn (T) sin —] cos -j— cos — 

x y 

where A mn (T) , B mn (T) are real and imaginary parts of slowly 
varying mode amplitudes; m and n are integers giving the number of 
half-lengths in the x— and y— directions; L x and L y are the sides of the 
experimental basin and wave vector k = (k x , k y ) = (2nm/L x , 2irn/L y ). 
First notice that for the case L x = L y = L the dispersion relation will 
take form 

u\k) = \gk+^ 

which means that taken above usual form of dispersion function 
corresponds to the unit square domain. The resonant solutions in a 
square with side L could be found as solution of resonance conditions 
with uj(k) written as a function of L. For instance, for \ijv = 16 in 
square basin with side L = 2 there is a resonant triad 

[1,14][23,13][24,27]; (25.0785, 50.2490, 75.3275). 

Notice that this triad is not resonant in the square basin with 
L = 1 and vice versa, the triad 
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[1,6] [4,5] [5,11]; (15.5681, 16.2945, 31.8626) 



is resonant for L = 1 and non-resonant for L = 2. Therefore, 
even in this very simple case when basin form does not change but 
only the basin sizes, no general selecting rule exists which allows us 
to find resonant triads in a given basin among the triads which are 
resonant in some other basin. Each time corresponding modification 
of dispersion function has to be used. For instance, for a rectangular 
basin the expression 

has to be used. 

Necessary experimental facilities 

A few wave generators must be available to make it possible to 
generate a few different wave frequencies simultaneously. In order to 
measure vertical deformation of the water surface from its quiescent 
position a wave gauge and a video recorder can be used. A wide variety 
of liquids does provide observable results. An experimental basin must 
have walls movable in such a way that, for instance, the basin form 
remains rectangular while the ratio of the basin sides changes or, in 
case of a circular basin, only the radius of the basin changes. 

4 Drift waves 

Till now our main interest was to investigate resonantly interacting 
waves because they play significant role in the energy transfer via the 
wave spectrum. But the opposite opposite situation - behavior of 
the waves which do not interact resonantly could also be of a great 
interest in some physical problems, for instance, for drift waves in 
Tokamak. The following questions arise: How long do these waves 
keep their energy and does the confinement time change from one 
mode to another? Does there exist a way to characterize the con- 
finement time for every specific wave? Does this time depend on the 
resonator form? etc. Evolution of drift waves in Tokamak plasma is 
described by Hasegawa-Mima equation which coincides up to renam- 
ing of variables with barotropic vorticity equation (BVE) which we 
used for numerical experiments described below. A possibility of en- 
ergy concentration in a given set of drift waves due to some instability 
mechanisms in plasma was conjectured by V.I.Petviashvili. Basing on 
this conjecture, we performed numeric simulations in which the most 
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part of the total initial energy of the wave field (till 60%) was located 
within a few specific wave groups while the rest of the energy was 
distributed randomly (Gaussian) or equally among all other waves of 
the spectrum. The results of numerical simulations show that all the 
waves can be partitioned into three classes having qualitatively differ- 
ent energetic behavior. 

Class A - active waves. It consists of the waves taking part in reso- 
nant interactions (RI) which are exact solutions of resonant conditions 
in general form 

u(k 1 )±uj(k2)±uj(k 3 ) =0, (4) 
fcl ± k 2 ± h = (5) 

and of the near-resonant waves which can be defined as follows. 
Let us denote 

R 1 = (Sj, Si Si), R 2 = (S 2 1: SlSi),...,RP = (S^SlS?) (6) 

the set of all wave vectors giving exact solutions of Sys. Q,© 
in some finite spectral domain T so that T contains all these wave 
vectors and some others which we denote as T r : 

T = {Sj}UT r Vi = l,2,..p, J = 1,2,3 

while r = 1, 2, 3, (T 2 — 3p ) where p < p depending on whether 
or not some of the triads R 3 have common waves. Let us say that 
waves are taking pari in ARI (approximate resonant interactions) if 
they have some non-zero frequency discrepancy 

< cj(fci) ±Lo(k 2 ) ±Lu(k 3 ) « (7) 

while Eq.© is hold exactly. Let us consider all waves taking part 
in ARI and choose the triads containing two waves of some specific 
resonant triad R J . Corresponding frequency discrepancy can be com- 
puted as 

w=± U (s})±u,(si)± U (r). (8) 

Waves TJ taking part in these interactions are called near-resonant 
waves. The class of active waves A is now defined as follows: 

A = {Sj}U7f CT, 
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The wave Tf o at which the minimum of IP is achieved for some 
specific j and vT r G T, is called a minimal near-resonant wave for 
the triad Energy leaves triad R? via the wave Tf . Using itera- 
tively this procedure of the construction of the minimal near-resonant 
wave for the triad (Sf, Sl,T r ) and so on, the way energy goes from 
one triad to another can be constructed. Energy oscillates first among 
the waves within resonant triad and then is re-distributed via a few 
near-resonant waves. The active waves transfer energy through the 
spectrum. 

Class P - passive waves. It consists of the waves taking part in ARI, 
except near-resonant waves. Therefore frequency discrepancy could be 
found as SV = uj(P t ) ± u{P k ) ± u{P r ) where Pj e T\{BP} and such 
Po does not exist for which some of waves (Pj&Pfc), (Pi&P r ), (Pfc&P r ) 
belong to resonant triad BP . Passive waves conserve their energy 
during time scale corresponding to the three-wave interactions and 
longer. The confinement time in some specific passive waves depends 
on the initial energy distribution of the wave field. For the same ini- 
tial conditions this time for different passive waves may vary 30-60% 
according to the magnitude of discrepancy. 

Class N - neutral waves. It consists of the waves which do not 
belong to class 1 or class 2, i.e. both Eqs.©,© are not satisfied, 
N = T\(A U P) and could be empty (see Table) depending on the 
form of resonator. But in case of N ^ the confinement time by 
these modes can be a few times longer than by passive modes. For 
illustration the graphs of the modes' energies are drawn with a relative 
shift by the axes Y. 

Numeric studies of BVE show different energetic behavior of the 
modes belonging to different classes. More precisely, energy of the 
neutral mode, is conserved during the time T conserve much longer then 
those predicted by wave turbulence theory, sometimes 3 to 5 times 
longer. Energy of the passive modes, with f2 = 0.21 and with 17 = 0.13 
conserves approximately during time T conserve /2 and afterwards begin 
to spread over the wave spectrum. Energy of the active mode, O = 0, 
is changing periodically from the very beginning. At the initial mo- 
ment T = all modes have the same amount of energy. 

Therefore, for given dispersion relation (i.e. for given resonator 
form and boundary conditions) in the finite spectral domain it is pos- 
sible to compute all its active, passive and neutral waves. The results 
for three different forms of resonators and two different spectral trun- 
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cations, T10 and T20, are following: 

• unit sphere 

— T10: Number of active modes is 4 and number of neutral 
modes is 3; 

- T20: Number of active modes is 51 and number of neutral 
modes is 3; 

• square 

— T10: Number of active modes is 15 and number of neutral 
modes is 0; 

— T20: Number of active modes is 53 and number of neutral 
modes is 0; 

• rectangular basin with side ratio 1/4 

- T10: Number of active modes is 4 and number of neutral 
modes is 75; 

- T20: Number of active modes is 16 and number of neutral 
modes is 300. 

Remark. The important fact is that the mode with given wave 
vector k = (m, n) may belong to different classes depending on the 
resonator geometry. For instance, the mode with wave vector (2,4) is 
active in a square and neutral in a rectangular basin with sides' ratio 
1/4. Therefore, if we supply energy to the mode with a given fre- 
quency, the time of energy confinement by this mode is defined by the 
resonator geometry and we can easily increase the number of passive 
waves or decrease the number of active waves just by changing the 
sides' ratio of the basin. The resonator geometry is indeed crucial fact 
in all these considerations because the neutral modes do not exist in 
some geometries and therefore the best confinement time principally 
can not be obtained in them. 

Our main hypothesis can now be formulated as follows 

The H-mode discharge in Tokamaks could possibly be described in 
terms of three-waves interactions of the drift waves. The H-modes are 
the neutral modes or passive modes with big frequency discrepancy, the 
L-modes are the passive modes with small frequency discrepancy and 
ELM-modes are the active modes. 

Let us construct a few parallels between the properties of H-mode 
discharge and the properties of the three-waves interactions in order 
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to support this hypothesis. We use the comprehensive review on H- 
mode discharges in Tokamaks (Hj to extract the most important facts 
on the subject. 

• The appearance of an H-mode discharge does not depend neither 
on the technique for its obtaining nor on the Tokamak geometry. 

- The three-wave interaction dynamics is defined completely by the 
initial energy distribution among the modes and therefore does not 
depend on way of "putting" the energy into the specific waves. The 
fact of this dynamics existence also does not depend on the Tokamak 
geometry. The crucial fact is the very existence of resonator which 
means the boundary conditions importance; as soon as modes have 
lengths long enough "to notice" boundary conditions the dynamics 
described above will take place. The specific geometry defines only 
which specific modes will be active, passive or neutral in this geometry. 

• The occurrence of repetitive instabilities at plasma edge - ELM-s. 
There exist ELM-s of the three types: "giant" (Type I), "grassy" 
(Type II) and Type III with big, small and intermediate frequen- 
cies correspondingly. 

- The active waves perform the periodic energy exchange between 
the modes of the resonant or near-resonant triads. Numerical simula- 
tions show that the periodic energy oscillations within the modes of 
these triads take place and their frequencies belong to the three differ- 
ent types described above. We do not touch here the question about 
specifics of energetic behavior of the passive waves in the presence 
of some active waves (it is in itself very interesting and complicated 
problem) and only point out the fact that under some conditions the 
presence of active waves may decrease the confinement time of some 
specific passive waves. 

• The H-mode is achieved only when a exceeded threshold in the 
heating power is exceeded but at the same time the improvement 
in core confinement is due to the reduction of turbulence. 

- These high and low boundaries for the total wave field energy 
form, so to say, "the corridor" of weak nonlinearity: the treshold in 
the heating provides the transfer from the linear regime to a weakly 
nonlinear one while reduction of turbulence prevents transition to the 
strong nonlinear regime. 

• The H-mode is initiated by the formation of the transport bar- 
rier; there exist a very strong temporal correlation between de- 
velopment of the transport barrier, a reduction in fluctuation of 
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amplitudes and a large change in the plasma edge; the mecha- 
nism by which the improvement in the core confinement is facil- 
itated by the transport barrier is not known. 

- The transport barrier formation leads to some specific profiles 
in the flow velocity and, therefore, to the specific energy distribution 
among the wave spectrum. Namely, the reduction of turbulence in the 
most part of the plasma core and increasing turbulence at the plasma 
periphery leads to the energy distribution needed for the resonant in- 
teractions become dominate, i.e. the most part of the wave field is 
settled into a few modes while the rest of energy is distributed more 
or less homogeneously among the others. If in an experiment we are 
able to increase the energy of some desired mode, we will obtain the 
theoretically predicted confinement. 

The presented analysis is neither full nor exhaustive and in any 
case has only instructive importance because the dispersion relations 
under consideration are not realistic for toroidal plasma. Our purpose 
here is just to attract attention to the complicated problem of finding 
realistic dispersion function(s) for this very interesting physical prob- 
lem. This hard common work of physicists and mathematicians may 
be well rewarded because as a result a new Tokamak geometry may 
evolve in which a lot of waves will improve their energy confinement 
time drastically. 



5 Concluding remarks 

We have shown that kinetic equation and clipping describe dynamics 
of two opposite limits of wave turbulence theory - continuous and dis- 
crete wave spectra correspondingly. 

The natural - and highly non-trivial - problem to be studied now 
is to describe transition from one limit to another. One of the most 
tricky point here is following. As it was noticed in Section 3, discrete 
systems described by the same evolution PDE, the same eigenfunc- 
tions of linear part of original PDE and having functionally the same 
dispersion relation will have different resonant triads corresponding to 
slightly different boundary conditions so that triad ([1,6] [4,5] [5,11]) is 
resonant in square basin with a side L = 1 (in non-dimensional units) 
and is not resonant if L = 2. 
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Barotropic vorticity equation demonstrates much more compli- 
cated behavior, namely even the form of linear modes itself changes 
when regarded with different boundary conditions. For instance, BVE 
on a sphere has a linear wave 

^sphere = AP™{sm.(j>) expi[mA + — — —rrt], (9) 

n{n + 1) 

where P™(x) is the associated Legendre function of degree n and order 
to; dispersion relation has then form 

^sphere = -2m/[n(n + 1)] 

and wave vector k = (to, n) has integer coordinates to, n. In the case 
of infinite plane BVE has a linear wave has form 

tppiane = A exp i(k x x + k y y + uj) (10) 

and wave vector k = (k x , k y ) has real coordinates k x , k y while disper- 
sion relation is written out as 

^plane = k x /(l + k x + ky). 

Thus, intuitively reasonable opinion, that it is enough to regard the 
resonance surface constructed for real variables and then find (if pos- 
sible) its integer points, does not hold here. Transition from a sphere 
to a special infinite plane called /3-plane can only be constructed lo- 
cally, at some interaction latitude and this transition is possible iff 
this latitude exists In particularly, it means that 1) transition 

from a spherical domain to an infinite plane is in fact transition to 
a one-parametric family of infinite planes; and 2) such a transition is 
not always possible. 

Nevertheless, these systems have much in common, not only gen- 
eral energetic behavior but also some more fine properties. For in- 
stance, "corridor" of weak nonlinearity described in previous Section 
and known from experimental plasma researches, has been established 

in numerical computations with capillary waves. A very promising 
idea has been formulated in - "to explain the slowdown of en- 
ergy flux and lower Kolmogorov constant observed by Pushkarev and 
Zakharov in numerical experiments |13j " in terms of discrete system 
behavior. It is a highly demanding task indeed to prove this hypothe- 
sis because transition from kinetic equation to the original PDE with 
a given boundary conditions will be much more complicated than just 
changing the interaction domain within the frame of the same equa- 
tion. 
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